Phase Transitions in the 1-d Long-Range Diluted Heisenberg Spin Glass 



O 

(N 



C/3 



I 

o 
o 



> 

(N 
rn 

rn 
O 



X 



Auditya Sharma^ and A. P. Young^'Q 

^Department of Physics, University of California, Santa Cruz, California 95064 

(Dated: January 15, 2013) 

We use Monte Carlo simulations to study the one-dimensional long-range diluted Heisenberg spin 
glass with interactions that fall as a power, a, of the distance. Varying the power is argued to be 
equivalent to varying the space dimension of a short-range model. We are therefore able to study 
both the mean-field and non-mean-field regimes. For one value of a, in the non-mean-field regime, 
we find evidence that the chiral glass transition temperature may be somewhat higher than the spin 
glass transition temperature. For the other values of a we see no evidence for this. 
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I. INTRODUCTION 

In the theory of phase transitions it is often helpful 
to study models in a range of dimensions ranging from 
above the "upper critical dimension" , du , where mean- 
field critical behavior is expected, to below the "lower 
critical dimension" , di , where fluctuations destroy the 
transition. However, it has been difficult to do this nu- 
merically for spin glasses, since d„ = 6 is quite large, and 
slow dynamics prevents more than a few thousand spins 
being equilibrated at low temperature T. It follows that, 
at and above cZ„, one cannot study a sufficient range of 
sizes to perform the necessary finite-size scaling analysis. 

As a result, there has been a lot of recent attention on 
long-range models in one-dimension, in which the inter- 
actions fall off with a power of the distance. Such models 
have a venerable history going back to Dyson [T,"?], who 
considered a ferromagnet with interactions falling off like 
and found a paramagnet-ferromagnet transition for 
1 < a < 2. Kotliar et al. \§\ studied the spin glass ver- 
sion of this model, which has received a lot of attention 
numerically in the last few years 

There are few analytical results on spin glasses beyond 
mean field theory. For the long-range models, Kotliar 
et al. [H computed critical exponents in an expansion 
away from the point where mean field theory occurs (e- 
expansion), but, as we shall see, this is poorly converged. 
Hence, most of what we know has come from numerical 
work. Much of this, including numerics on long-range 
spin glass models, has studied the Ising version. However, 
there are also reasons to study models with vector spins, 
such as the Heisenberg (3-component) model. 

One motivation is that Kawamura @ proposed that 
there are two separate transitions in vector spin glasses, 
a spin glass transition at T = Tsg and a "chiral glass 
transition" at higher temperature Tec involving a freez- 
ing of vortex-like variables called chiralities. While the 
original scenario had Tsg — 0, it now appears that Tsg 
is non-zero in three or more dimensions, but the ques- 



tion of whether Tsg < Tqg, or whether there is a single 
transition at which both types of ordering occur, is still 
open 

A second motivation for studying the Heisenberg, 
rather than Ising, spin glass is that it is possible to study 
larger sizes, see for example Ref. [19j, which should be 
helpful in a finite-size scaling analysis. In a second pa- 
per in this series [2l|, we will investigate whether there 
is a de Almeida- Thouless [121 (AT) line of transitions in 
a magnetic field for Heisenberg spin glasses. This follows 
our recent work psj which shows that there is an AT line 
for vector spin glasses provided one considers a random 
field. The ability to study larger sizes will be particularly 
useful for the AT-line study. 

Here we present data for the zero field transition for 
the Heisenberg spin glass for values of the parameter a 
corresponding to (i) the mean-field regime, (ii) the non- 
mean-field regime, and (iii) the borderline case where the 
transition disappears. Most of our results find no evi- 
dence for separate spin-glass and chiral-glass transitions. 
However, for one set of parameters in the non-mean-field 
regime, the data indicates that Tqg niay be somewhat 
greater than Tsg- Whether this result remains valid in 
the thermodynamic limit, will require future studies on 
significantly larger systems. 

The plan of this paper is as follows. In Sec. |ll] we de- 
scribe the model that we study, while in Sec. lIIII we gives 
some technical details of the simulations. The results are 
presented in Sec. llVl and our conclusions summarized in 
Sec. El 



II. MODEL 



We consider the Hamiltonian 
H - 



(1) 
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where Si 
length 1, 



are classical 3-componcnt Heisenberg spins of 
and the interactions are independent vari- 
ables with zero mean and a variance which falls off with 
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a power of the distance between the spins, 
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ij\av 
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where [• • • ]at) means an average over disorder. In the ver- 
sion used in early studies Q], every spin interacts with 
every other spin with a strength which falls off, on av- 
erage, like Eq. ([2]). However, this means that the CPU 
time per sweep varies as N"^ , rather than iV, so large sizes 
can not be studied. This problem was solved by Leuzzi 
et. al. ^ who proposed a model in which, instead of 
the magnitude of the interaction falling off with distance 
like Eq. ([5]), it is the probability of there being a non-zero 
interaction between sites which falls off, and when 
an interaction does exist, its variance is independent of 
rij. The mean number of non-zero interactions from a 
site, which we call z, can be fixed, and we take z = 6. To 
generate the set of pairs that have an interaction 

with the desired probability we choose spin i randomly, 
and then choose j (7^ i) at distance r.y with probability 



P^j 



(3) 



where, for rij, we put the sites on a circle and use the 
distance of the chord, i.e. 



N . 
— sm 

TT 



(4) 



If i and j are already connected, we repeat the process 
until we find a pair which has not been connected before. 
We then connect i and j with an interaction picked from a 
Gaussian interaction whose mean is zero and whose stan- 
dard deviation is J, which set equal to 1. This process is 
repeated precisely Nb = zN/2 times. 

The result is that each pair {i,j) will be connected 
with a probability which must satisfy the condition 
NJ2j Pij = since Pij only depends on |i — j|, Pa — 0, 
and there are precisely Nz/2 connected pairs. It follows 
that, for a fixed site i, 



[ Jij \ av — J Pij — J 



(5) 



Note that Pij is different from z xpij in Eq. ^ because of 
the constraint that no bond can occur twice. The mean- 
field spin glass transition temperature for m-component 
vector spins is given by [23 | 



1/2 



rriMF 
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where the last equality follows from Eq. ([5]) . We set J = I 
so that, for the situation here. 



J = l, 



6, m = 3 , 



(7) 



TABLE L A summary of the behavior for different ranges of a 
in one space dimension and at zero field. Infinite range means 
that {j^i) Jij diverges unless the bonds Jij are scaled by 
an inverse power of the system size. The behavior is mean- 
field like for a < where — 2/3, and a finite-temperature 
transition no longer occurs for a > ai where cr; = 1. 



a 




Behavior 


(7 = 




SK model 


< 0- 


< 1/2 


Infinite range 


1/2 < 


a < 2/3 


Mean-field with Tsg > 


2/3 < 


a < 1 


Non-mean-field with Tsg > 


1 < a 




Tsg = 



we have 



rpMF _ 

^SG — 
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0.816. 



(8) 



the same as for the nearest-neighbor Heisenberg spin 
glass on a simple cubic lattice. 

By varying a one finds different types of behav- 
ior d, [13, as shown in Table HI For < 1/2 the model 
is non-extensive (for instance the mean-field transition 
temperature in Eq. ^ diverges) unless the interactions 
are scaled by an inverse power of the system size. We 
will call this "infinite range" . The extreme limit of this, 
CT = 0, is the Sherrington-Kirkpatrick i26i] model, whose 
exact solution was found by Parisi |27h29| . In fact, it has 
been suggested [soj (see also Ref. [sil) that, in the ther- 
modynamic limit, the behavior of the model is identical 
to that of the SK model for the whole range < a < 1/2. 

The model is extensive for cr > 1/2 and a finite temper- 
ature transition is expected for cr < cr/, where the "lower 
critical" value is 



ai = l. 



(9) 



The transition is in the mean-field universality class Q 
for cr < cr„, where the "upper critical" value is 

au = 2/3 . (10) 

For au < CT < o'l, there is a finite-temperature transition 
with non mean-field critical exponents. In this paper we 
will study both mean-field and non mean-field regions. 
Finally for a > ai the there is no transition at finite 
temperature. 



III. NUMERICAL SETUP 

We perform large scale Monte-Carlo simulations for 
cr = 0.6,0.75,0.85 and 1. From the previous section we 
note that cr — 0.60 is in the mean-field regime, cr = 0.75 
and 0.85 are in the non mean-field regime, and cr = 1 
is the borderline case, a = ai, beyond which there is 
no transition. A plausible scenario is that Tsg = for 
cr = 1, though the possibility that Tsg is non-zero cannot 
be ruled out a priori. Table HIl lists the parameters of the 
simulation. 
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TABLE II: Parameters of the simulations. A'samp is the num- 
ber of samples, A^cquii is the number of overrelaxation Monte 
Carlo sweeps for equilibration for each of the 2Nt replicas 
for a single sample. The same number of sweeps is done in 
the measurement phase, with a measurement performed every 
four overrelaxation sweeps. The number of heatbath sweeps 
is equal to 10% of the number of overrelaxation sweeps. Tmin 
and Tmax are the lowest and highest temperatures simulated, 
and Nt is the number of temperatures used in the parallel 
tempering. 



cr 


N 


3amp 


TV -1 




7m ax 


Nt 


0.6 


128 


16000 


128 


0.20 


0.70 


40 


0.6 


256 


16000 


256 


0.20 


0.70 


40 


0.6 


512 


16000 


512 


0.20 


0.70 


40 


0.6 


1024 


16000 


1024 


0.20 


0.70 


40 


0.6 


2048 


16000 


2048 


0.20 


0.70 


40 


0.6 


4096 


6100 


4096 


0.30 


0.70 


40 


0.6 


8192 


1000 


8192 


0.30 


0.70 


50 


0.6 


16384 


500 


16384 


0.35 


0.70 


55 


0.6 


32768 


400 


32768 


0.35 


0.70 


60 


0.75 


128 


8000 


128 


0.20 


0.55 


40 


0.75 


256 


8000 


256 


0.20 


0.55 


40 


0.75 


512 


8000 


512 


0.20 


0.55 


40 


0.75 


1024 


3000 


1024 


0.20 


0.55 


40 


0.75 


2048 


3000 


2048 


0.20 


0.55 


40 


0.75 


4096 


3000 


4096 


0.20 


0.55 


40 


0.75 


8192 


1100 


8192 


0.20 


0.55 


50 


0.75 


16384 


500 


16384 


0.25 


0.55 


55 


0.75 


32768 


400 


32768 


0.25 


0.55 


60 


0.85 


128 


16000 


512 


0.09 


0.30 


40 


0.85 


256 


16000 


1024 


0.09 


0.30 


40 


0.85 


512 


10000 


2048 


0.09 


0.30 


40 


0.85 


1024 


10000 


8192 


0.14 


0.22 


20 


0.85 


2048 


8000 


16384 


0.14 


0.22 


20 


0.85 


4096 


4000 


32768 


0.14 


0.22 


36 


0.85 


8192 


2000 


65536 


0.15 


0.21 


20 


0.85 


16384 


1700 


131072 


0.16 


0.21 


20 


1.0 


128 


2000 


2048 


0.03 


0.10 


10 


1.0 


256 


2000 


4096 


0.03 


0.10 


10 


1.0 


512 


2000 


16384 


0.02 


0.10 


40 


1.0 


1024 


1200 


524288 


0.017 


0.08 


40 


1.0 


2048 


500 


2097152 


0.017 


0.08 


60 



A. Equilibration 

As discussed in earlier work [5','32l] there is a convenient 
test for equilibration with Gaussian interactions, namely 
the relationship 



J2 2 



Is) , 



(11) 



is valid in equilibrium but the two sides approach 
their common equilibrium value from opposite di- 
rections as equilibrium is approached. Here U — 
-(l/-^)E(i,i) '^ijJiji^i ■ Sj)]av is the average energy per 
spin, qi = {l/Nb) J2{i,j) • Sj)^]a^v is the "link over- 

lap", and qs = (l/^^b) [((Si ' Sj)2)]av, where 

Nb = zN/2, and ey = 1 if there is a bond between i and j 
and is zero otherwise. Equation is easily derived by 



integrating by parts the expression for the average energy 
with respect to Jy since it has a Gaussian distribution. 
Note that in the numerics we set J = 1. 

We determine both sides of Eq. pT|) for different num- 
bers of Monte Carlo sweeps (MCS) which increase in a 
logarithmic manner, each value being twice the previ- 
ous one. In all cases we average over the last half of 
the sweeps. We consider the data to be equilibrated, if, 
when averaging over a large number of samples, Eq. (jlip 
is satisfied for at least the last two points. 



B. Simulation Technology 

To equilibrate the system in as small a number of 
sweeps as possible, with the minimum amount of CPU 
time, we perform three types of Monte Carlo sweeps [TtI - 

M- 

The workhouse of our simulation is the "Microcanon- 
ical" sweep [S^ (also known as an "over-relaxation" 
sweep). We sweep sequentially through the lattice, and, 
at each site, compute the local field on the spin, — 
J'ijSj. The new value for the spin on site i is taken 
to be its old value reflected about H, i.e. 



S' 



(12) 



These sweeps are microcanonical because they preserve 
energy. They are very fast because the operations are 
simple and no random numbers are needed. For reasons 
that are not fully understood, it also seems that they 
"stir up" the spin configuration very efficiently 18| and 
the system equilibrates faster than if one only uses "heat- 
bath" updates, described next, see e.g. Fig. 9 of Ref. js^ . 

We also need to do heatbath sweeps in order to change 
the energy. As for the microcanonical case, we sweep se- 
quentially through the lattice. We take the direction of 
the local field H^, to be the polar axis for the spin on 
site i. We compute the polar and azimuthal angle of the 
new spin direction relative to the local field by the re- 
quirement that this direction occurs with the Boltzmann 
probability, see Ref. for details. 

Finally we perform parallel tempering sweeps [ ssl . l36j 
to prevent the system from being trapped in local min- 
ima at low temperature. We take Nt copies of the system 
with the same bonds but at a range of different tempera- 
tures. The minimum temperature, T,„iii = Ti, is the low 
temperature where one wants to investigate the system 
(below TsG in our case), and the maximum, Tmax = Tjv^ , 
is high enough that the the system equilibrates very fast 
(well above Tsg in our case) . A parallel tempering sweep 
consists of swapping the temperatures of the spin config- 
urations at a pair of neighboring temperatures, Ti and 
Ti_|_i, for i = 1,2,- •• jTnt-i with a probability that 
satisfies the detailed balance condition. The Metropo- 
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lis probability for this is [3£ 



We define the chiral correlation length by 



P(r swap) 



exp(A/3A£;), (if A/3 A£; < 0), 
1, (otherwise), 



(13) 

where A/3 = l/T.-l/T.+i and AE = Ei-Ei+i, in which 
Ei is the energy of the copy at temperature T,;. In this 
way, a given set of spins (i.e. a copy) performs a random 
walk in temperature space. 

We perform one parallel tempering sweep for every ten 
overrelaxation sweeps. Since there are two copies of spins 
at each temperature, indicated by labels "(1)" and "(2)" 
in Eq. (|14p below, we actually perform parallel tempering 
sweeps among the set of Nt copies labeled "(1)" and, 
separately, among the set of Nt copies labeled "(2)". 



C. Quantities Measured 

The main quantities measured in this simulation are 
the spin glass susceptibility XSG, and the chiral glass 
susceptibility XCG, at wavevectors k — 0, and k = 2n/N, 
and from these we obtain the two corresponding correla- 
tion lengths, ^sG and ^cG- The spin glass order param- 
eter, q^'^{k), at wave vector k, is defined to be 



1 ^ 



m(i) 



•^(2) ik-Ri 



(14) 



i=l 



where fi and v are spin components, and "(1)" and "(2)" 
denote two identical copies of the system with the same 
interactions. We run two copies of the system at each 
temperature in order to evaluate quantities such as the 
spin glass susceptibility, defined in Eq. dTSl) below, with- 
out bias. From this we determine the wave vector depen- 
dent spin glass susceptibility xsGik) by 



XsG{k)=Nj2m'"'ik)\')]s 



(15) 



where (• • • ) denotes a thermal average and [•••]av de- 
notes an average over disorder. The spin glass correlation 
length is then determined from 



^SG 



1 / Xsg{0) 

2sin(fcniin/2) V XSG (^min) 



1/(2^-1) 



, (16) 



where fcmin — (27r/L). For the Heisenberg spin glass, 
Kawamura defines the local chirality in terms of three 
spins on a line as follows [lo| : 

Ki = S,;+i • Si X Si_i. (17) 

The chiral glass susceptibility is then given by 

XCGik) ^ N[{\q,{k)\^)U, (18) 
where the chiral overlap (7c (k) is given by 



N 



(19) 



GG 



Xcg(O) 
2sin(fc,nin/2) \ 



- 1 



l/(2<7-l) 



(20) 



As will be revealed in the next section, three of the four 
quantities defined above, XsGi£,sGi and XCG may be used 
in a finite-size-scaling analysis to locate and analyze the 
phase transition. 



D. Finite-Size Scaling 

According to finite-size scaling [s^l, the correlation 
length of the finite-system varies, near the transition tem- 
perature Tc, as 



N 



= X[N^'%T - T,)] , (2/3 < fT < 1), (21a) 
= X[N^'^{T - Tc)] , (1/2 < cr < 2/3), (21b) 



in which the correlation length exponent, is given, in 
the mean-field regime, \yy v = 1/(2ct — 1). We will use 
Eq. (PT|) for both the spin glass correlation length ^sgi 
in which Tc will be set to Tsg, and the chiral glass cor- 
relation length S^cGi in which will be set to Tea- It 
follows that, if there is a transition at T = Tc, data for 
^/iV {^/N"^^ in the mean-field region) for different sys- 
tem sizes N should cross at Tc. 

We also present data for xsG = Xsg(O), which has the 
finite-size scaling form 



XSG 



= C[N^/-{T-T,)] , (2/3<a< 1), 



(22a) 



^ C[N'/^{T ~ Tc)] , (1/2 <a< 2/3). (22b) 

Hence curves of Xsg/^^ (xsg/^^^"^ in the mean-field 
regime) should also intersect. This is particularly use- 
ful for long-range models since r\ is given by the sim- 
ple expression 2 — 77 = 2(t — 1 exactly. However, we do 
not know the exponent corresponding to 77 for the chiral 
glass susceptibility, so we will not use this quantity in the 
finite-size scaling analysis. 

In practice, there are corrections to this finite-size- 
scaling, so data for different sizes do not all intersect 
at the exactly the same temperature. Including lead- 
ing corrections to scaling, the intersection temperature 
T*(iV, 2iV) for sizes N and 2Af varies as 



A 



(23) 



where A is the amplitude of the leading correction, and 
the exponent A is given by 



A = - + w 

V 



(24) 



where uj is the leading correction to scaling exponent. 



5 
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FIG. 1: (Color online) The main figure shows data for 
^sa/N^^^ , in which the power of A*' is chosen following 
Eq. (|21bp with v = l/(2cr — 1), as a function of T for different 
system sizes for a = 0.6. The inset shows data for xsc/N^^^ , 
in which the power of A'^ is chosen following Eq. (I22bl) . 



FIG. 3: (Color online) The main figure is a log-log plot of the 
logarithmic derivative of £^sG for a = 0.6 for different sizes A'' 
evaluated at the intersection temperatures T* (TV, 2N) shown 
in Fig. [2] According to Eq. (|21b|l the slope is expected to be 
1/3. The best fit is a little smaller than this, indicating that 
corrections to scaling are still present for these sizes. The 
inset is the same but for XSG- 
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FIG. 2: (Color online) The intersection temperatures, 
T*{N,2N), obtained from the data in Fig. [H for ^sg/N^^^ 
and Xsg/N^^'^ for a = 0.6, as a function of A'^ ^, with 
A — 0.467 determined from Eq. (|24p (which is valid in 
the MF regime). A fit using all 8 points from ^sG gives 
Tsg = 0.564 ±0.002, while a fit using the data for the largest 
5 pairs of sizes from xsG gives Tsg = 0.562 ± 0.002. 



IV. RESULTS 

A. a = 0.6 (mean-field regime) 

As shown in Table U cr = 0.6 lies well inside the mean- 
field regime. Hence, according to Eq. (|21bl) . results for 
£,Sg/N^/'^ should intersect at Tsg with v set equal to 



1/(2(7 — 1). The data is shown in the main part of 
Fig. [TJ The intersections do not occur at precisely the 
same temperature, but fitting the intersection temper- 
atures to Eq. is helped by the fact that we know 
A = I — 2(7 in the MF regime [4l|, which gives a value 
0.467 here. A straight line fit of £,sg/N''/'^ against N'^, 
shown in Fig. H gives Tsg = 0.564 ± 0.002. 

The inset to Fig. [1] shows data for Xsg/N^^^, which 
should also intersect at Tsg according to Eq. (j22b[) . This 
time, we find that corrections to scaling are well described 
Eq. ([^51) but only if we consider just the largest five pairs 
of sizes. The fit to Eq. shown in Fig.H gives Tsg = 
0.562±0.002, which is consistent with that obtained from 
the spin glass correlation length. 

We have also measured the chiral glass correlation 
length. However, we find that ^cg/N"/^ < 10"^^ in 
the vicinity of Tsg- Hence chiralities can not play an 
important role in this range of temperature. 

According to Eqs. (|21b[) and (|22bp the argument of 
the scaling functions is N^/^{T — Tc). Hence, at Tsg 
the, logarithmic derivative of S^sG a-nd XSG should vary 
as N^/^ . As we have seen, the data do not all intersect at 
the same temperature, and so we evaluate the derivatives 
at the intersection temperatures T* plotted in Fig.O The 
results are shown in Fig. [S] We get a power of 0.28 from 
^SG and 0.29 from xsG, in both cases a little less than 
1/3, indicating that there are still some corrections to 
scaling even for these large sizes. 
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FIG. 4: (Color online) The main figure shows data for differ- 
ent sizes for £,sg/N at a — 0.75 which, according to Eq. (|21a|) . 
should intersect at Tsg- The inset shows data for xsg/N"^^^ 
(with 2 — ?7 = 2a — 1), which should also intersect at Tsg 
according to Eq. (|22ap . 




We fit the intersection temperatures to Eq. ([23]) . but 
unfortunately we do not know the value of the exponent 
A outside the MF region, and have to include it as fit 
parameter. The fit is therefore to a non-linear function 
of the parameters. We determine the fit p arameters using 
the Levenberg-Marquardt algorithm [42| . The data from 
XSG is better behaved than the data from ^sG so we use 
the former to determine the exponent A and then fix this 
value in the fit (which did not include the smallest size) 
to the data from ^so- This procedure is justified since 
the exponent giving the leading correction to scaling, A, 
is universal, though the amplitude of this correction {A 
in Eq. ((23)) ) is non-universal. The results are A = 0.44 ± 
0.13, Tsg = 0.359 ± 0.003 from ^sg, and Tsg = 0.354 ± 
0.005 from xsG- The two estimates for Tsg agree within 
the error bars. 

The data for £_cg/N in the region of the spin glass tran- 
sition temperature is very small, around 10~*. Hence, as 
for a = 0.6, chiralities do not play an important role in 
the transition. 

According to Eqs. (|21al) and (|22ap . adapted to include 
corrections to finite-size scaling, the logarithmic deriva- 
tive of CsG and xsG should vary as N^/" at T*{N, 2N). 
The plots in Fig. [5] yield l/vsG = 0.25 from ^SG and 
0.29 from XSG- The difference presumably comes from 
corrections to scaling. 

Kotliar et al. J| calculated critical exponents to leading 
order in e = tr — 2/3 with the result 



^-3(1 



12e + 0(e2)) . 



(25) 



The large coefficient of e indicates the expansion becomes 
poorly converged well before the "lower critical" value 
e = 1/3 (fT = 1). Even for the present value of tr, which 
corresponds to e = 1/12, Eq. ([25]) gives 1/v — Q, and so 
is not useful for comparison with the numerics. 



C. (J = 0.85 (non mean-field regime) 



FIG. 5: (Color online) The intersection temperatures for cr = 
0.75. The all 8 data points from xsG were fitted to Eq. (|23|) . 
with the result that A = 0.44 ± 0.13. The same exponent was 
then used to fit the results for the largest 7 (pairs of) sizes 
from CsG, for which sub-leading corrections to scaling seem 
to be more important. The resulting values for the transition 
temperature are Tsg = 0.354 ± 0.005 from xsG, and Tsg ~ 
0.359 ± 0.003 from Csc 



B. a — 0.75 (non mean-field regime) 

For a — 0.75 we are no longer in the MF regime. 
Hence, according to Eqs. (|21a|) and (I22ap . ^sg/N and 
Xsg/N'^~^ (with 2 — 77 = 2(7—1), should intersect at 
Tsg- The data are shown in Fig. U and the resulting 
intersection temperatures are shown in Fig. [Sj 



For a = 0.85 we are further in the non-mean-field re- 
gion. Among the different models studied here, this is 
the one which is most similar to a short-range model in 
three dimensions. The spin glass data is shown in Fig. [T] 
Chiral correlations are larger than for a — 0.6 and 0.75, 
so we show data for S^cG in Fig. IS] 

As for a = 0.75, to extrapolate the intersection tem- 
peratures to the thermodynamic limit, we resort to 
Levenberg-Marquardt fits with three parameters. Using 
data from xsG we find Xsg = 0.99 ± 0.13. The data 
from £^cG is insufficient to determine Acg since the fits 
give Acg — 0.79 ± 0.74, i.e. the error bar is as large as 
the best estimate. Nonetheless, in estimating Tcg and 
its error bar from the data for xcGj we allow A to vary. 
For there appear to be very large sub-leading correc- 
tions, so we fixed the value of Xsg to that obtained from 
XSG when fitting the results from S,SG- 
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FIG. 6: (Color online) The main figure is a log-log plot of the 
logarithmic derivative of ^sG for a = 0.75 for different sizes N 
evaluated at the intersection temperatures T*{N, 2N) shown 
in Fig. [S] According to Eq. (|21al) the slope is expected to be 
1/v. The inset is the same but for XSG- 



FIG. 8; (Color online) Data for £,cg/N, the spin glass cor- 
relation length divided by system size, as a function of T for 
different system sizes for a = 0.85. 




0.15 0.16 0.17 0.18 0.19 0.2 0.21 

T 



FIG. 7: (Color online) The main figure shows data for ^sg/N, 
as a function of T for different system sizes for a = 0.85. The 
inset shows data for xsg/N^~'^ with 2 — = 2cr — 1. 



Results for r*(A;', 2N) against and fits are shown 

in Fig. ini In the plot, for all data we use the value of 
A determined from xsG- However, we again emphasize 
that, in the fit to the XCG data, we estimated Tcg and 
its error bar by allowing A to vary. From the fits, we find 
TsG = 0.167±0.001from xsa, Tsa = 0.166 ±0.004 from 
tsG, and Tcg = 0.190 ± 0.006 from ^co- 

The two estimates of Tsg agree with each other but 



are lower than Tcg- This would imply spin-chirality de- 
coupling for a = 0.85. However, we note that the data for 
the spin glass correlation length appears, at intermediate 
sizes, to be extrapolating to a value for Tsc of around 
0.19 (our value for Tec) but then, for the largest sizes, 
veers down to about 0.167 (very close to our value for 
Tsg obtained from xsg)- Hence we cannot rule out the 
possibility that a similar "crossover" may occur for the 
chiral glass correlation length data, but at even larger 
sizes. If so, then spin-chirality decoupling would not oc- 
cur at the largest scales. We also note that the actual 
values of £,cg/N shown in Fig. |5] are still very small in 
the vicinity of Tcg, about 1/50 of the value of £,sg/N 
around the transition, see Fig. [T] Hence we are very far 
the regime with ^cG > S,SG which will ultimately occur 
in the presence of spin-chirality decoupling. 

Figure [TO] shows results for the logarithmic derivative 
of ^SG and xsG evaluated at T*{N,2N). Fits give 1/z/ 
= 0.22 from £^sG and 0.29 from xsG- The curvature in 
the data for ^SG indicates strong finite-size corrections. 
Presumably these corrections are also the reason for the 
difference between the two estimates for l/iy. Figure 
[TT] shows similar data but for £,cg- The best fit gives 
l/v — 0.21, which is close to the estimate from £,sg- 
Note that the coefficient of N^^'^, 14.01, is much larger 
than the corresponding value, 3.18, for ^sg, presumably 
to compensate for the overall value of ^cG being much 
less than that of ^sG in the vicinity of the intersection 
temperatures T*. 
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FIG. 9: (Color online) The intersection temperatures for a — 
0.85. The data from XSG was fitted to Eq. (|23p . with the result 
that A — 0.99±0.13. This is the value used to plot all the data. 
The same exponent was then used to fit the largest four sizes 
in the data from ^SG, for which sub- leading corrections appear 
to be very large. The full Levenberg-Marquardt fit to the 
data for ^ca was unable to determine A with any precision, 
see text. However, in estimating Tea and its error bar from 
data for ^cG we allowed A to vary. The resulting values for 
the transition temperatures are Tsa = 0.167 ± 0.001 from 
XSG (all 7 points), Tsa = 0.166 ± 0.004 from ^sa (4 points), 
and Tea = 0.190 ±0.006 from ^cg (5 points). We emphasize 
that, in these estimates, we fixed the value of A only for the 
data from ^sG- 



25 



20 - 



15 - 



10 



1 1 — 1 — 1 — 1 

a = 0.85 


1 1 1 1 1 


1 1 


— 1 — 1 — r- 


12' 












3.18 n"^^ - 




,x 
















y'A 


200 


1000 


10000 






. 1 


1 






16 






,A-' 




N)±| 




,A- 






1- 

:2 8 




X 




A 


C3 
W 
?^ 










O 


'A 


1.21 n''-2=' 1 




? 4 






ml 






, ,.l 



100 



1000 



10000 



FIG. 10: (Color online) The main figure is a log-log plot of the 
logarithmic derivative of ^sg for a = 0.85 for different sizes A'' 
evaluated at the intersection temperatures T*{N, 2N) shown 
in Fig. [§1 According to Eq. (|21al) the slope is expected to be 
1/v. The inset is the same but for xsg- 



90 
80 

i 70 
V 

60 - 
'g 50 - 

O 

V 40 - 



30 



- 

— 






1"\ 

a = 0.85 








A''' _ 


- 


















---''a 












14.0 N°-^' 


A 






..1 1 



100 



1000 
N 



10000 



FIG. 11: (Color online) Similar to Fig.[TO]but for ^cG- 
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FIG. 12: (Color online) Data for ^sg/N, the spin glass cor- 
relation length divided by system size, as a function of T for 
different system sizes, for a — 1.0. The inset shows the in- 
tersection temperatures T*{N, 2N), as well as a log- log fit as- 
suming TsG = 0. The fit works well for iV = 128, 256 and 512. 
There is no intersection for the two lowest sizes, A'' — 1024 
and 2048 for T greater than the lowest temperature we could 
simulate, 0.017. This temperature is well below the value of 
the fit extrapolated to N = 1024, which is about 0.042. Hence 
the intersection temperatures actually fall off faster at large 
sizes than shown in the fit. 



D. 



It is known from the early work of Kotliar et al. Q 
that a = 1 is the "lower critical" value <ti, above which 
there is no spin glass transition. Interestingly, Viet and 
Kawamura |20| claim that chiral glass ordering persists 
to slightly larger values of a. Testing this claim is one 
of our main motivations for performing simulations at 
cr = 1. 

Figure [T^ shows the finite scaling for ^SG- In the inset. 
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FIG. 13: (Color online) Data for £,cg/N, the chiral glass cor- 
relation length divided by system size, as a function of T for 
different system sizes, for a = 1.0. The inset is a log-log fit to 
the intersection temperatures T*{N, 2N), assuming Tea = 0. 
The data is quite consistent with this behavior. 
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FIG. 14: (Color online) Data for xsg/N, the spin glass sus- 
ceptibility divided by system size (the power 2 — 77 is equal to 1 
here), as a function of T for different system sizes for a — 1.0. 
There are no intersections for the temperature-range in which 
the simulation is conducted. 

we show a log- log plot of T*(A^, 2iV) versus N and include 
a straight-line fit for N = 128, 256 and 512. This fit works 
well. We find no intersections in the range of T that we 
can equilibrate (T > 0.17) for N = 1024 and 2048. Hence 
the data is well consistent with Tsg — 0. 



Figure [13] shows the corresponding figure for £,cG- 
Again, the inset shows log-log fits assuming that the tran- 
sition temperature. Tea in this case, is zero. The fit is 
satisfactory and indicates that we do not find a finite 
value for Tcg at ct = 1, in contrast to the claim of Viet 
and Kawamura (20l |. 

Figure [HI shows the data for XSG- There are no inter- 
sections at all in the range of temperature that we can 
equilibrate, consistent with the conclusion from the ^so 
data that Tsg = 0. 



V. CONCLUSIONS 

Our primary motivation to study the Heisenberg spin 
glass in one dimension with long-range interactions which 
fall off with the power of the distance, is to test Kawa- 
mura's spin-chirality decoupling scenario in which Tea > 
TsG: and his subsequent prediction [13] that chiral glass 
ordering persists for a > ai, where tr; = 1 is the "lower 
critical" value for the spin glass transition, with a finite 
value of Tcg at a — I. 

For (T = 1 we find Tcg — Tsg = in contrast to Viet 
and Kawamura [20]. For most of the other values of a 
our data is well consistent with a single phase transition. 

However, for a — 0.85 the best fits for the sizes we 
can equilibrate indicate Tcg > Tsc, see Fig. [9} Inter- 
estingly that figure shows very strong subleading correc- 
tions to finite-size scaling in the data for ^SG since it 
only fits Eq. ([23]) for the largest sizes. At intermediate 
sizes the intersection temperatures seem to be heading 
towards the chiral glass transition temperature obtained 
from the ^SG data, but then dip down for the largest 
sizes. We therefore cannot rule out that similar behavior 
might occur for the chiral glass correlation length but at 
even larger length scales. In this case there would be no 
spin-chirality decoupling. We also note that, for a given 
size and temperature, ^cG remains considerably smaller 
than ^SG in the vicinity of the intersection temperatures 
T*, compare the main part of Fig. [7] with Fig. [5] Hence 
the data is very far from the regime with £^cG > isG 
which ultimately prevails for Tsg < T < Tec if there is 
spin-chirality decoupling. 

In a subsequent paper [2l| , we will investigate, for the 
same models, under what circumstances an AT line of 
transitions occurs in a magnetic field. 
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